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o 

■ The BPS Skyrme model is a specific subclass of Skyrme-type field theories which 

O I possesses both a BPS bound and infinitely many soliton solutions (skyrmions) satu- 



rating that bound, a property that makes the model a very convenient first approx- 
imation to the study of some properties of nuclei and hadrons. A related property, 
the existence of a large group of symmetry transformations, allows for solutions of 
rather general shapes, among which some of them will be relevant to the description 
p I of physical nuclei. 

pH I We study here the classical symmetries of the BPS Skyrme model, applying them to 

construct soliton solutions with some prescribed shapes, what constitutes a further 
. important step for the reliable application of the model to strong interaction physics. 

o\: 

^ ■ PACS numbers: ll.30.Pb, 11.27.+d 
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INTRODUCTION 



The Skyrme model [l| (SM), a non-linear field theory for an S'f/(2)-valued field, is meant 
to be a low energy effective theory, describing some interesting aspects of strong interaction 
I physics. In this model, pions play the role of primary fields (excitations around the trivial 

vacuum), whereas nucleons and nuclei are, on the other hand, represented by topological 
solitons, collective excitations which are part of the nonperturbative spectrum of the theory. 

The application of the SM to nuclear and hadronic physics has been quite successful at a 
qualitative level Q-Q, but it encounters some problems once a more detailed, quantitative 
agreement, is required. The main obstacle for this is the absence of (almost) BPS solutions 
in the original SM, as well as in its standard generalizations. Indeed, although there exists 
a BPS bound already in the original model, as proposed by Skyrme, nontrivial soliton 
solutions cannot saturate this bound. As a consequence, higher solitons, meant to describe 
larger nuclei, are strongly bound, in striking contrast to the weak binding energies of physical 
nuclei. 

Some alternatives approaches to improve this situation have been recently advanced. 
Basically they imply the extension of the symmetries of the Skyrme type theory to conformal 
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transformations [6| or to volume preserving diffeomorphisms 7|. It it the aim of this paper 
to further elaborate on one of them, namely, the proposal of 7|. 

The SM may be generalized in a rather straightforward way, by simply adding some judi- 
ciously chosen extra terms to its defining Lagrangian j^- 12|. Indeed, the addition of extra 
terms becomes a quite natural step when one recalls the fact that the SM is an effective 
theory, supplemented with the condiment of some simplicity and symmetry constraints. In 
fact, assuming, as one usually does, that we want to maintain the field content of the orig- 
inal model, as well as its Poincare invariance and the standard Hamiltonian interpretation 
(Lagrangian quadratic in time derivatives), the number of possible terms is in fact quite re- 
stricted. One may then just have: a potential term (no derivatives), a standard kinetic term 
(the nonlinear sigma model term) quadratic in first derivatives, the 'Skyrme term' originally 
introduced by Skyrme (quartic in derivatives) and, finally, a term which is the square of the 
baryon number current (topological current), which is sextic in derivatives. 

As it has been demonstrated in there is a submodel, termed 'BPS Skyrme model' 
(BPSSM), defined by a Lagrangian consisting of just the potential and sextic terms, which 
satisfies some quite interesting properties. Indeed, it possesses a BPS bound, and infinitely 
many BPS solutions saturating this bound. Besides, it has been also shown in 0] that the 
static energy functional of the model is invariant under an infinite number of symmetry 
transformations, a fact that is obviously related to the properties enunciated in the previous 
sentence. 

Among the symmetry transformations, an interesting type are the volume preserving 
diffeomorphisms (VPDs), since they are precisely the symmetries of an incompressible fluid, 
a fact pointing to a possible relation to the liquid drop model for nuclei. The BPSSM, 
therefore, has several appealing features from the point of view of the description of nuclei 
(see, for example, UM)- The model is, in fact, constructed assuming that the coherent 
(topological) excitations play an especially important role in strong interaction physics. This 
assumption is directly related to the suppression of the usual kinetic term in the Lagrangian 
and, as a consequence, one might expect that the BPSSM will not lead to reliable results 
in the weak field regime. To overcome these shortcomings, it may be necessary to augment 
the lagrangian by further structures for a more consistent description of nuclear or hadron 
physics. There are, for instance, initial data for which the BPSSM does not have a well- 
defined Cauchy problem; thus, a standard kinetic term must either be added explicitly or 
induced by quantum corrections, to remedy this situation. 

We think that the BPSSM provides an approximation which may be quite reliable for 
the study of static properties and for the dynamics in a region of relatively high density 
(i.e., with a not too small baryon charge density) like, e.g., in a soliton background. On the 
other hand, it will generally not be reliable in near-vacuum regions, and moreover cannot be 
applied at all to consider perturbative phenomena corresponding to quantum fluctuations 
of the pion field around the vacuum, since the dominant term would then be non-quadratic. 

Because of the above, it would be important to relate the properties of solutions of the 
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BPSSM to the corresponding solutions of more general Skyrme-type models. A stumbling 
block which immediately pops up when attempting this task is the different sizes of the 
respective spaces of solutions, which are in turn due to the different symmetry groups of 
the models. The solutions of the BPSSM may have almost any symmetry, due to the 
huge symmetry group of the field equations. In particular, there are spherically symmetric 
solutions (i.e., with spherically symmetric energy densities) for all the possible values of the 
baryon charge, Qb- This is not the case, on the other hand, for the original SM and its 
non BPS generalizations. Typically, the = 1 skyrmion is spherically symmetric, the 
Qb = '2 skyrmion has cylindrical symmetry, while higher-charge skyrmions have, at most, a 
set of discrete symmetries. Indeed, their energy densities are invariant under some discrete 
subgroup of the rotation group SO (3) (see, for example, [l3]-[l6|). A skyrmion of the BPSSM 
with the same set of discrete symmetries would, therefore, be a good starting point for the 
inclusion of physical effects induced by adding other extra terms to the Lagrangian. Because 
of this, it would be important to find a method for the systematic construction of solutions 
of the BPS Skyrme model with some prescribed symmetries. 

It is the purpose of the present notes to investigate the space of BPS solutions further, 
making explicit use of its symmetries as a tool to generate new solutions. To that end, we 
shall take the spherically symmetric ones as a starting point for the construction. Finally, 
we shall show that all local solutions may, in fact, be constructed in this way. 

This article is organized as follows: In section [Til we define the model, and introduce our 
notation and conventions. Then we construct the classical Hamiltonian in section lllli The 
BPS bound is considered in IIVI In |V] we explore the issue of symmetries, within both the 
Lagrangian and Hamiltonian contexts. In |Vl] we derive and discuss the main properties of 
the BPS solutions of the model. We also construct several explicit classes of solutions with 
some prescribed symmetries, including the important case of discrete symmetries. In IVIII 
we summarize our results and conclusions. 



II. THE MODEL 



The Lagrangian density £, which has an SU{2) valued field U as dynamical variable, 
may be written as follows: 

£ = £o6 = -A^^ ^^^M _ ^2^^u^ f/t) ^ (1) 
where A is a positive constant, B^^ denotes the topological current: 

B'' = ^ e'"""'tT{UL,L^) , = U^d.U , (2) 

and V is a potential density. The current 5^ is 'topologically conserved', namely, it can 
be shown to be conserved, regardless of the equations of motion. The resulting conserved 
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charge, Q_b, is therefore given by: 

Qb = J d'xBo = d^xe'Hi{LiL^Lk) 

= ^ j d'xe''Hi{U^dm^djUU^dkU) , (3) 

the degree of the map — )■ S^, an integer which is invariant under arbitrary globally 
well-defined coordinate transformations, as well as under global isospin rotations of U. It 
is, in fact, invariant under the much bigger group of target space transformations leaving 
invariant a certain target space volume form, see below. 

To proceed to the classical equations of motion, it is convenient to introduce a specific 
parametrization for the three degrees of freedom of U. 

Following 0] , we use a real scalar field ^ plus a 3-component unit vector fi, so that: 

U{x) = e^«(^')"(^)-" , (4) 

where r are the three Pauli matrices. The real scalar ^ runs from to vr, while the two 
independent parameters defining n may be taken as the two components of a complex 
variable u, by means of a stereographic projection: 



n 



(u + u, —i{u — u) , \u\'^ — l) . (5) 



i + \u\ 

In this way, one obtains for the Lagrangian density an expression in terms of C,, u and u: 

^ = TTxrw ^^^p^^f - -"'^(0 (6) 

(1 + |m| j 

where the lower indices in those variables denote partial derivatives with respect to the 
spatial coordinates, and we have assumed that the potential may only depend on U through 
tiU. 

With the notation = d^V, the Euler-Lagrange equations read: 
where 

These objects satisfy, by construction, the relations 

H^u^' = H^u^ = 0, = K,n^ = 0, H^^^ = K,u^ = l^e^^^'^ i^u^u^f , (8) 



fyr^^ ) = , (7) 



which are often useful. 
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III. HAMILTONIAN AND STATIC ENERGY 



In order to construct the Hamiltonian, we first introduce a more compact notation, in 
terms of tliree real fields with a = 1, 2, 3, such that u = ^^^^ + i^^'^^ and ^^^^ = ^. 
Then C may be written as follows {^^""^ = do^'^"'''): 

where the kinetic term is determined by a metric Gi^ab)-, given by: 

(-^(ab) - (^(1))2 + (^(2))2]4 y^^) 

where: 

Qt^ = e.,ke'^''if\t\ (11) 

In order to see whether the system defined by C is regular or not, we note that Q = [Q^f^] 
the 3x3 matrix defined by the nine elements Q-"^ {i = 1, 2, 3; a = 1, 2, 3) is proportional to 
the cof actor matrix of the matrix X = [■Cl"'']- 

Q = 2 cof (X) . (12) 

Thus, we see that the metric [G(afe)] (hence, the Lagrangian system) is regular if and only if 
det[^f"^] 7^ 0. In other words, the regularity of the system is equivalent to the non vanishing 
of the Jacobian determinant: 

J = det[X] = det [^] 7^ , (13) 

OXi 

for the mapping between the sphere (i.e., one-point compactified M^) in coordinate space 
and the one in SU{2). 

Under the assumption that ( |T3l) holds true, the inverse of G = [G'(afe)] may be found by 
elementary algebra. Indeed, 

IC^H^b) _ [1 + (^(^)^ + (^^^^)T ^(a)^(&) (.,^ 
^ J 8A2j2sin4(^(3)) • 

Thus, the Hamiltonian density in terms of the variables ^^'^\ its spatial derivatives, and 
their canonical momenta n^'^-', becomes: 

^ [1 + (^W)^ + (,) 
16A2 J2 sin4(^(3)) 
4 \2 j 4 f (3) ( ijkAl)A2)A3)\ 2 
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which, for a Lagrangian hke the one we are considering, coincides with the energy den- 
sity of the system. In particular, for the static configuration case to be considered in the 
forthcoming sections, the total energy E is: 

We have shown that the regularity of the system depends on the field configurations 
considered. Specifically, the system is singular in regions where the fields take their vacuum 
values (.^'^'^^ = const, such that V(.^''^'*) = 0). This already demonstrates that, while the 
system may provide a good approximation to the description of static properties of nucleons 
and nuclei via solitons (Skyrmions) and for the dynamics in regions with nonzero baryon 
charge density (where it is regular by construction), its fully consistent apphcation to dy- 
namical nuclear physics requires additional structures like, e.g., quantum corrections, or the 
inclusion of further terms in the Lagrangian. 



IV. BPS BOUND 



The static energy functional in f|T6|) . or, in terms of the variables and u introduced 
previously, 

' A^sin^e 



E = j d^x 
obeys a Bogomolny bound. Indeed 
Asin^e 



'i + \u\^Y 



(17) 



E 



d'^x 



mni-^ - , A7^^ f ,3 2/iAsin^^A/V 

e <mMnM/ ± /iV Vj T J d^^X ^^^^^2^2 <mUnUl 



\u 



2\2 



+ u 



\u 



2\2 



± (2A/i7r^ 



-i /■ 3 sin^^A/V 

Ti 1 + m 



2A/i7r2(VV)|5| 



(18) 



where {\/V) is the average value of \/V on the target space S'^. The corresponding Bogomolny 
(first order) equation is 



Asin^e 



(19) 



(1 + |n|2)2 

The static second order field equations may be derived from the squared Bogomolny equation 
by applying a gradient dk and by projecting onto eijkdj^^^^ where ^^"'^ = u, u). We remark 
that a completely analogous BPS bound can be found for the BPS baby Skyrme model in 
one lower dimension 17j- 20|]. 
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Another interesting observation is that the BPS equation can be formulated in the lan- 
guage of a non-linear generalization of the static (vacuum) Nambu-Poisson equation. Indeed 



the left hand side can be recast into the Nambu-Poisson three-bracket 21 1 



{^..,«.,C}^.™.,«1^«! (20) 

where the target space embedding coordinates X^^ A = 1, 2, 3, 4 form a three-sphere (i.e., 
(X^)^ = 1) and are related to the previous coordinates like X" = n^sin^, a = 1,2,3, and 
X^ = cosC,- Then, the generalized Nambu-Poisson dynamics is given by 

^ = e^sc-D (x^,X^,X^} +X^v/V(X^, (21) 

(JjL 



which differs from the standard case by the additional factor \/V in the last term 21 |. 
Obviously, although the dynamics of the BPS Skyrme model is profoundly different, the 
BPS equation provides static solutions to this generalized Nambu-Poisson equation. Such 
solutions may be interpreted as vacuum configurations of the underlying hyper-membrane 
Lagrangian [22|. We remark that if one assumes from the outset that the target space 
variables X^ span a three-sphere, as we do in this paper, then there is no dynamics in Eq. 
fl^ . i.e., = 0, as follows from the fact that the r.h.s. of ( 12T1) is proportional to X"^ in this 
case. This just corresponds to the well-known result that the static vacuum equations for 
;he hyper-membrane imply that the brane embedding coordinates X"^ span a three-sphere 



22l |. So, our model generalizes the static hyper-membrane action, with a correspondence 
between the BPS solitons and the vacuum membrane configurations, but with completely 
different dynamics. 

It may be instructive to compare the BPS bound arising above with a 1 + 1 dimensional 
analogue: the search for (non-trivial) static minimum energy configurations for the Sine- 
Gordon model. Here, a real scalar field (p is in the presence of a potential density U{(f) 
which allows for non-trivial topology. The Lagrangian density is: 

C = \{d,vf-U{^) (22) 



nr r . , a/a 



,4 



U{v) = [l-cos(— ^)] . (23) 
A m 

The static energy is then: 

/ + 00 T 
dx,[-{d,ify + Uiv)] . (24) 

The non-negative potential has non-trivial minima for 

277" 777 

^ = ^^ = -^N, iVeZ, (25) 
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all of them having zero energy. Finite energy vacuum configurations must tend to one of 
the minima when Xi — ±00. 

The topologically conserved current is j'^ = ^^^^"duV i.l^^''^ = 1)2), which obviously 
satisfies 9 ■ j = 0. Its associated topological charge is quantized: 

Q^ = ^ dx.dM^) = N , (26) 

it is a constant of motion, and it is akin to a winding number, if one interprets (f as an 
angular variable. 

Note the striking similarity with the BPS Skyrme model, when one writes the energy as 
follows: 

E, = / dx,[U'-^njof+Ui^)] . (27) 



-2 A 

The static energy then also verifies a Bogomolny-like bound, since: 



(28) 



Thus: 



E,>±V2 r^dx.^^vw) 

J —oo UXi 



3 

2v^7r^|g^|. (29) 
A 



where: 



m\ = [ (30) 



the average of \/lA{}p) over the fundamental region. 

Of course, the first order equations that result from saturating the bound may be found 
by other methods; they lead to the well-known static solutions by a single quadrature. What 
we learn from the comparison with this model is that the particular form of the Lagrangian 
of the BPS Skyrme model involving the square of the topological current, is what makes 
it produce quite powerful constraints on the solution. It is interesting to note that the 
kinetic term in this 1+1 dimensional example allows for two different interpretations, either 
as a standard kinetic term or as the topological current squared, which is no longer true 
in higher dimensions. In other words, the simple Sine-Gordon type soliton model in 1+1 
dimensions allows for two different generalizations to higher dimensions, generalizing either 
the standard kinetic term or the topological current, and the model studied in the present 
paper just corresponds to the second case. 
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V. SYMMETRIES 

The Lagrangian certainly has the standard Poincare symmetries. Besides, the sextic term 
is the square of the pull back of the target space volume form on S^, 

dV = —i . ^" I ^.^.^ d^dudu (31) 

(1 + \u\^Y 

so this sextic term is invariant under target space diffeos which do not change this form 
(the volume preserving diffeos (VPDs) on S'^). The potential only depends on ^, so it is still 
invariant under those diffeomorphisms which do not change ^, i.e., under the diffeos which 
obey 

^ — ^ , u^u{u,u,^), (1 + ImI^)"^^^^^^^ = (1 + ImI^)"^^^^^^. 

The symmetries mentioned so far are symmetries of the action, i.e. Noether symmetries. 

The static energy functional has some further symmetries. Indeed, it is invariant under 
volume preserving diffeos on the base space M^, as can be seen easily. The Bogomolny 
equation has even more symmetries as we want to demonstrate now. For this purpose we 
introduce the new target space coordinates 

u = ^e'* = tan(x/2)e^^ H{g) ^ 



(for later convenience we also introduced which together with ^ and $ provides the 
standard hyperspherical coordinates on the target S^), and 

and rewrite the Bogomolny equation as 

VF(0 ■ VH{g) X V$ = ±1 (33) 
or, in terms of differential forms 

dFdHd(f) = ±dx^dx^dx^ (34) 

from which it is obvious that the Bogomolny equation has as its symmetries all the VPDs 
both in base space and in a modified target space defined by the volume form dFdHd^. The 
above equation implies, in fact, that all local VPDs on base space produce local solutions 
of the BPS equation. The problem is that, in general, a local solution cannot be extended 
to a global one, because of the different geometry and topology of the base space and the 
modified target space. The modified target space is defined by the volume form 

dFdHd^ = ; — ^ sin xd^dxd^ (35) 
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and differs from the volume form on S'^ by tlie additional factor \/\/V. There does not 
exist a unique riemannian metric giving rise to this volume form, but a natural choice which 
assumes that the spanned by u (i.e., x ^"^^ ^) remains intact is 

ds' = de + ^?^(^^' + si^' Xd^')- (36) 

For V = 1 this is just the round metric on in hyperspherical coordinates, but for nontrivial 
potentials the resulting target space manifold is different. Indeed, potentials which may 
support finite energy skyrmion solutions must have vacua = where V(^o) = 0, and 
the above metric is singular at the vacuum values ^o- These singularities may either be 
integrable (i.e., the function F defined in (132|) is well-defined and finite even at vacuum 
values ^ = ^o); in which case the total volume of the modified target space is still finite. 
In the opposite case, the total volume is infinite. One further conclusion may be drawn 
immediately by integrating Eq. ( I34l) . If the total volume of the modified target space is 
finite, then any skyrmion solution of the BPS equation must have compact support (i.e., be 
a " compacton" ) . Further, its volume must be equal to \B\ times the total volume of the 
modified target space, where B is the winding number. For equivalent results for the case 
of the BPS baby Skyrme model in one lower dimension, we refer to 20|. 
We remark that for V = sin^ ^ the metric on the target space describes in fact a 3 dimensional 
cylinder with a very simple skyrmion solution (see below). 



VI. SOLUTIONS 

As already said, locally, any VPD on base space will provide a solution of the BPS 
equation, but this solution will, in general, not be extendible to a global, genuine one (i.e., a 
skyrmion), because of the nontrivial topology one should have on the modified target space. 
A more promising strategy is the following: start from a simple known solution which may 
follow from a simple ansatz. Then one may generate new solutions by composing the given 
solution with a VPD on base space M?. If the VPD is well-defined on the whole of M^, 
then it will map genuine skyrmions into genuine skyrmions. In the case of compactons, we 
may even relax this condition, since it is then sufficient for the VPD on base space to be 
well-defined in the region of the compacton. 

To proceed, let us first find some simple solutions with the help of an ansatz in spherical 
polar coordinates 

e = e(r), X = X{0), ^ = n^ (37) 
which inserted into the BPS equation yields: 

A sin^ ^ 

sin dxd^ = =i=r ^ sin 9drd9dip , (38) 
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leading to x = ^ and 

_r;_sin^^^ = Tr^dr 
or, after the coordinate transformation: 



y = -^^r^ (39) 



to the autonomous ODE: 

= (40) 



y2V(!) 

We have chosen the sign which leads to a negative ^y, which is compatible with the boundary 
conditions ^(r = 0) = vr, ^(r = oo) = for a potential which takes its vacuum at = 0. 

Let us consider now the symmetries of these solutions. This issue depends on the criterion 
used to characterize that symmetry. Note that a given solution will not be invariant under 
any rotation, because it depends on the two angular coordinates 6 and tp. The energy density, 
on the other hand, depends only on the radial coordinate r and is, therefore, spherically 
symmetric. Note, however, that there exists another symmetry criterion, often used for 
solitons, whereby there is spherical symmetry when the effect of a base space rotation on a 
solution can be undone by a corresponding target space rotation. Under this criterion, only 
the solution with topological charge n = 1 is spherically symmetric (i.e., all rotations can 
be undone). Solutions with higher winding number n only have cylindrical symmetry, i.e., 
only a rotation about the z axis ip (p + a can be undone by a target space rotation (a 
phase transformation u — e"*""^). 

In any case, we shall call all solutions of the spherically symmetric ansatz "spherically 
symmetric solutions" in what follows. We shall first review some general properties of 
these spherically symmetric solutions and, in a next step, construct solutions with lesser 
symmetries. 



A. Solutions with spherical symmetry 

Many qualitative aspects of solutions maybe easily derived from the particular form of the 
potential, which should be contrasted with the typical situation in general Skyrme models, 
where similar results usually require a full three-dimensional numerical simulation. 
First of all, depending on the form of the potential in the vicinity of the vacuum, one can dis- 
tinguish three types of solitonic configurations: compactons (where the solution approaches 
its vacuum value at a strictly finite distance) and exponentially as well as power-like local- 
ized solutions. Using the BPS equation and expanding the potential at a vacuum (e.g., at 
^ = 0), V = Vo^" + one easily finds that for a < 6 one gets compactons. There is also 
one exponentially localized solution for a = 6, while for a > 6 we find power-like localized 
solitons. 



12 



Another important feature of solutions reflects the number of vacua of the potential. It is 
easy to prove that for one-vacuum potentials the BPS solutions are of the nucleus type (no 
empty regions in the interior), while two- vacuum potentials lead to shell-like configurations. 
Let us present some particular examples. For the most elaborated family of one vacuum 
potentials, the so-called old potentials 

V,i,= (tt(^^]Y ^ V(0 = (l-coser (41) 



(where a is a real positive parameter), we find (besides the previously known compacton) a 
solution with exponential tail (a = 3) in implicit form 

cos ^ + 111 tan ^ ~ ~ ^ 
, and power-like localized solutions. E.g., for a = 6 we get 



^ = 2 arc cot Y 3V2y. 
A family of two-vacuum potentials is given by 

V..e./=(Tr(l^)Tr(i±^)y ^ 1/(0 = (1 - cos^ ^ (42) 

which is the chiral counterpart of the so-called new baby potential. The vacua exactly 
coincide with the boundary values for the scalar field i.e., = 0, vr. From the BPS property 
of the solution one can immediately see that the energy density should have a shell structure 
with two zeros: one at the center of the soliton, while the second (outer zero) can be located 
at a finite distance (compact shells) or approached asymptotically at infinity. Without losing 
generality (the potential is symmetric under the change of the vacua) we assume that ^ = 
is the outer vacuum. Of course, the inner vacuum can only be reached at a finite point as 
y > 0. This implies that only compact solitonic shells are acceptable. Specific examples of 
exact solutions are, for a = 1 



and for a = 2 



arccos(y2y - 1) y e [O, V2] 

y>V2, 



n — 



IT - V2y y G 

y>T2- 

The latter solution is, in fact, a solution for the case when the target space is a three- 
dimensional cylinder, as sin^^/^A^ = const. 

In order to deal with non-compact shell skyrmions, we need to modify our potential in such a 
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way that one vacuum (say, the inner vacuum at ^ = vr) is always approached in a compacton 
manner . A simple choice is 



hell II 



Tr 



l + U 



Tr 



1 - U 



y(0 = (l + cosO(l-coser (43) 



1- (2^ 



Again, we find compact shell skyrmions a < 3 

{arccos 


an exponentially localized skyrmion for a = 3 

^ = = arccos 



3-a . 
V2 ' 



y<& 



z > 



V2_ 
3— a 



1 - 2e" 



-V2y 



and shell skyrmions which extend to infinity but are localized in a power-like manner (a > 3) 



^ = arccos 



3-a a - 3 ^ 
2— + ^^y 



V2 



B. Solutions with cylindrical symmetry 

Now we assume that a spherically symmetric solution has been found, and we want to 
use symmetry transformations to map them to new solutions. In a first step we construct 
solutions with cylindrical symmetry, using the ansatz (in cylindrical coordinates) 

^ = ^{p,z), g = g{p,z), <^ = nif (44) 

where = (x^)^ + (x^)^, z = x^. The Bogomolny equation for this ansatz may be written 
hke 

dF^'^UH = ±dqdp (45) 

where F'-"^ = nF and 



9 = P = z 



or like the Poisson bracket 



{F^-\H} ^ - = ±1- (46) 

oq oq op oq 



Further, we know that it has the spherically symmetric solution 



g = g, = tan(e/2) = _^ = = 9s{q.p) (47) 

p'^ + z'^ + z ^JZq + p^+p 

and (depending on the potential) 

e = Ur) = UV^Q + ^) = Uq^p)- (48) 
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As a consequence, a general solution with spherical symmetry may be written like 

='^.(Q(g,p),^(g,p)), 9iq,P) = 9siQiq,p),P{q,p)) (49) 

where {Q,P) are related to {q,p) via a canonical transformation, i.e., {Q,P} = 1. 
A first class of examples is given by 

Q = U{q), P ^ 



U'{q) 



where U'{q) 7^ V g must hold. Further, it should hold that liiRg^Q U (q) / q = const, to 
have a well-behaved function near p = 0. Among these examples the scale transformation 
Q = a'^q, P = a^^p can be found, which corresponds to the scale transformation — >■ ax^, 
x'^ — >■ ax"^ and x^ — )■ a~'^x^. Another class of examples is 

C. Solutions with discrete symmetries 

Here, we want to construct a class of base space VPDs which transform solutions with 
spherical or cylindrical symmetry into solutions which only preserve symmetries w.r.t. to 
some discrete rotations about the z axis. Concretely, we want to consider solutions which 
may be written like 

^ = ^{p^z) = ^s{p,z), g = g{p,z) = gs{p,z), ^ = n(f (50) 

where $,s, gs, ^ = nip constitute a known solution with either spherical or cylindrical sym- 
metry. That is to say, we consider base space VPDs which act nontrivially only on p and ip, 
where for simplicity we restrict ourselves to the following transformations, 

P = p{p,'f), ^ = ^{^)- (51) 

Using q = as before, and q = q{q, (p), the condition for the transformation to be a VPD 
simplifies to 

dqdip = dqd(p. (52) 

A class of formal solutions is given by 

^ = fi¥>) (53) 



in close analogy to the results of the last section. In order to define genuine diffeomorphisms, 
however, the transformations have to obey some further conditions. In particular, for the 
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new coordinates q and ip to define polar coordinates on they must satisfy tlie boundary 
conditions 

g(g = 0,v2) = 0, g(g = cx), V?) = oo, 
^{ip = 0) = 0, (f{(f = 27r) = 271. (54) 

In addition, tlie vector field generating the flow induced by the coordinate transformation 
must be well-defined (nonzero and nonsingular) on the whole of M^. A class of examples 
fulfilling all the required conditions is given hj f = (f + (c/m) sinm</), i.e., by the class of 
transformations 

g = (1 + ccos'mip)~'^ q m G N 

ip = (p -\ sinmo9 c G M, IcI < 1. (55) 

m 

Clearly, if a solution ^i"'\p, z, if) is invariant under rotations about the z axis (in the sense 
that its energy density is invariant under these rotations), then the new solution ^i"'\p, z, (p) 
is invariant only under the discrete set of rotations + (27r/m). 



VII. SUMMARY 

We explored in detail the symmetries of the static energy functional of the BPSSM, and of 
its related BPS equation. Then we applied these symmetries to the systematic construction 
of new solutions, start ing from known ones. This is in the spirit of the dressing methods of 



classical integrability 23|], which is an open problem for higher dimensional generalizations 
24| . an initial motiviation of this work. Specifically, this allowed us to construct solutions 
with some prescribed symmetries, what is quite relevant to the physical problem one wants 
to consider. We gave concrete examples of solutions with cylindrical symmetry and with 
symmetries w.r.t. some discrete subgroup of the group 5*0(2) of rotations about the z axis. 
In this context, it would be interesting to construct solutions with the symmetries of platonic 
bodies or other discrete subgroups of the full rotation group 5*0(3) (crystallographic groups), 
because solitons with these symmetries frequently show up as true minimizers of the energy 



in the original Skyrme model or some of its generalizations ij]- 16|] . The corresponding 
volume-preserving diffeomorphisms producing solutions with these symmetries will be more 
complicated than the ones constructed in the present paper, and it almost certainly will be 
more difficult to find them. 
This issue is under current investigation. 
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